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A recent lattice study of the Landau-gauge overlap quark propagator has shown a connection 
between centre vortices and dynamical chiral symmetry breaking in SU(3) gauge theory. We fur¬ 
ther investigate this relationship through an exploration of the connection to the instanton degrees 
of freedom. After identifying centre vortices on the lattice in maximal centre gauge, we smooth 
configurations using multiple algorithms. We are able to create an instanton liquid-like background 
on configurations consisting solely of centre vortices, analogous to that found on Monte-carlo gen¬ 
erated configurations after similar smoothing. Through calculations of the static quark potential 
and Landau-gauge overlap propagator, we show that this background is able to reproduce all salient 
long-range features of the original configurations. Thus we conclude that the information necessary 
to recreate the long-range structure of SU(3) gauge theory is contained within the centre vortex 
degrees of freedom. 


I. INTRODUCTION 

It has long been known that a centre vortex m model 
of the QCD vacuum is capable of explaining confinement, 
and recent work [9] hints that this is the model of con¬ 
finement most consistent with lattice results. In SU(2) 
lattice gauge theories, centre vortex gauge configurations 
have been found to recover the full string tension, and 
thus confinement [TOl [H] . Centre vortices have also been 
found to be responsible for dynamical chiral symmetry 
breaking UMB, and thus are the dominant long-range 
structure of the SU(2) vacuum. 

In SU(3) gauge theory, however, fewer studies exist, 
and so the picture is less clear. While removal of centre 
vortices still removes the string tension - and thus con¬ 
finement - a background consisting solely of centre vor¬ 
tices reproduces just 66% of the original string tension 
[53] . Additionally, while studies of dynamical chiral 
symmetry breaking through the Landau-gauge AsqTad 
quark propagator were unable to show a loss of dynamical 
chiral symmetry breaking coincident with vortex removal 
[24| . this did occur in a study of the low-lying hadron 
spectrum . Recently, the work of Ref. , studying 
the Landau-gauge quark propagator using the chirally 
sensitive overlap fermion action, revealed removal of dy¬ 
namical chiral symmetry breaking following vortex re¬ 
moval, and furthermore reproduction of dynamical mass 
generation after small amounts of cooling on configura¬ 
tions consisting solely of centre vortices. 

These results raise the possibilty that previous stud¬ 
ies of centre vortex configurations in SU(3) gauge theory 
suffered from a lack of smoothness; since gauge field con¬ 
figurations consisting solely of centre vortices have links 
consisting entirely of centre elements of SU(3) by defini¬ 
tion, they are high action, rough configurations. In this 
work we will examine the effects of smoothing vortex- 
only configurations, comparing to smoothed Monte Carlo 


generated (“untouched”) configurations. We find a re¬ 
production of the instanton-like degrees of freedom on 
vortex-only configurations, providing a mechanism for 
the dynamical mass generation seen in Ref. |26j . 

We begin in section |T^ with a brief overview of the 
process of identifying centre vortices on the lattice. In 
section m we examine the instanton content of these 
configurations through direct examination of the gauge 
field. We perform the smoothing procedures of cooling 
and over-improved stout link smearing on the configura¬ 
tions, and examine the resulting instanton-like objects. 
Instantons are associated with low-lying modes of the 
Dirac operator EH, and thus through the Banks-Casher 
relation |28j provide a mechanism for dynamical chiral 
symmetry breaking. We show an essentially equivalent 
background of instanton-like objects emerging on both 
vortex-only and untouched configurations through the 
process of smoothing. In section jlVj we analyse heavy 
quark confinement using the static quark potential, and 
find, again, almost perfect agreement between vortex- 
only and untouched configurations after smoothing. In 
section we analyse dynamical mass generation, and 
thus dynamical chiral symmetry breaking, through the 
Landau-gauge overlap quark propagator, showing that 
vortex-only configurations are able to reproduce dynam¬ 
ical mass generation at all levels of smoothing considered, 
over a range of bare quark masses. 

II. IDENTIFYING CENTRE VORTICES ON 
THE LATTICE 

The strategy for identifying thick centre vortices on 
the lattice is simple; we seek to decompose gauge links 
in the form 

Uf,{x) = Zf,{x)R^{x), ( 1 ) 

where Z^{x) is an element of the centre of SU(3), 
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Z^{x) = exp 


( 2 ) 
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in such a way that all vortex information is captured in 
the centre-projected elements Z^{x), and R^{x) thus con¬ 
taining the remaining short-range noise. The centre flux 
through plaquettes P^i,{x) is identified via their centre 
projection, 

Pfiu{x) = Z^{x) Z„{x + fi) Zl{x + v) Zl{x), (3) 

where there exists a non-trivial centre flux if P^,v{x) ^ I. 
The set of such plaquettes form closed surfaces; the 
“thin” centre vortices. These are embedded within the 
thick centre vortices of the original, Monte Carlo gener¬ 
ated, configurations. 

The natural choice m to perform this decomposition 
is to find some gauge transformation Q{x) which min¬ 
imizes the distance between links and centre elements; 
i.e.. 
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'^WU^ix) - Zf,{x)\\ min, (4) 

noting that we minimize over both the gauge transforma¬ 
tion V,{x) and choice of centre element Z^{x). U^{x) de¬ 
notes the link U^{x) after performing the gauge transfor¬ 
mation ^l{x). Choices of gauge fixing functional fulfilling 
this condition are collectively known as Maximal Cen¬ 
tre Gauge. Since all centre elements are mapped to the 
identity in the adjoint representation, this is equivalent to 
maximising the trace of the links in the adjoint represen¬ 
tation. There are many possible choices [ini El m m- 
|32] of gauge fixing functional which achieve this. We use 
the so-called “mesonic” maximal centre gauge [32] . 

^ dim ^ 

X,fJ, 

where the normalisation factor, 1/ is intro¬ 

duced to guarantee |i?|< 1. Setting = I everywhere 
except at a single lattice site x', we locally maximize the 
quantity 

i?iocai(x') =EM|Trr!(x')C/^(x')P+ (6) 

j:,\TTU,{x'-f^)n\xT. 

We restrict il.{x') to be in one of the SU(2) subgroups of 
SU(3). We parametrize an SU(2) matrix [f2(a;')]su(2) by 

[f^(a;')]su(2) = 1 - o-j, (7) 

and embed the resulting SU(2) matrix in one of the 3 
SU(2) subgroups of SU(3). Then, i?iocai(a:0 can be re¬ 
written as 

4 ^ 4 

77local{4^ ) — ^ ^ 2 ^ ^ “b C. (8) 

i,j=l i 

Where are elements of a real, symmetric matrix, bi a 
real vector and c a real constant, all dependent only on 
Ufj,{x') and U^{x' — fi). A list of these coefficients can 


FIG. 1. Values of the maximal centre gauge fixing functional, 
R, after gauge fixing on a configuration-by-configuration ba¬ 
sis. Square points in blue provide values from the original 
ensemble generated in the Markov chain. Circular points in 
red correspond to values found after performing a random 
gauge transformation before gauge fixing. 


be found in the appendix. The global maximum of this 
quantity can be found using, e.g., the method of Ref. [33]. 
We then update the relevant links with the SU(3) gauge 
transformation ^l{x'), and iterate over lattice sites. Each 
iteration over SU(2) subgroups and then lattice sites con¬ 
stitutes a single sweep. We perform 20,000 such sweeps 
on each configuration, after which the maximum relative 
change in i?iocai(a;) over a single sweep was found to be 
less than 10“^°. 

We note that although we have globally maximised 
Riocai(a^) for n(x), we have only found a local maximum 
of R. We illustrate this Gribov copy issue in Fig. 
we show the values of R obtained on the same config¬ 
urations fixed to maximal centre gauge starting from 
two distinct random gauges. Although we have reached 
similar values of R both times, there is no correspon¬ 
dence on a configuration-by-configuration basis; variation 
is well outside of the stability of the gauge fixing algo¬ 
rithm, indicating we have reached distinct local maxima 
on each occasion. This issue means we are unable to 
uniquely identify the vortex matter of a configuration; 
some of the physical centre vortices embedded in the lat¬ 
tice may be missed by the vortex identification procedure. 
Herein, we have selected the original ensemble generated 
in the Markov chain, denoted by the square blue points 
in Fig. [2 

While choices of gauge fixing functional exist which 
avoid the issue of Gribov copies, i.e. Laplacian centre 
gauge |131[3D], the vortices found in this way do not scale 
appropriately in the continuum limit, and so lack an in¬ 
terpretation as physical objects [22]. By contrast, the 
vortex density of objects found through Mesonic Maxi¬ 
mal Gentre Gauge fixing have a density independent of 
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lattice spacing, and thus admit a continuum interpreta¬ 
tion [221 [24]. 

After gauge fixing, we decompose links into polar form 
as 

TiUj^{x) = rf,{x) exp[ie^{x)], (9) 

and thus select the centre element Z^{x) = 

exp[^m^(x)], m^(x) G {—1,0,1}, such that 
is closest to 0fj_{x). 

Results are calculated on 50 pure gauge-field configura¬ 
tions using the Liischer-Weisz 0(a'^) mean-field improved 
action [34] , with a 20^ x 40 volume at a lattice spacing of 
0.125 fm. In Fig. we show a histogram of the values of 
d^{x) found after fixing to maximal centre gauge, show¬ 
ing three clear peaks, corresponding to the three centre 
phases. 



0 

FIG. 2. Histogram of the phases of the trace of links on a 
sample configuration before and after fixing to maximal centre 
gauge. Values after MCG fixing are denoted by the red circles, 
and random gauge by blue squares. 

We thus define the three ensembles which will serve as 
the basis for this work: 

• the untouched configurations, 

C/^(x); (10) 

• the vortex-only configurations, 

Z^{x); ( 11 ) 

• and the vortex-removed configurations, 

R4x) = Zlix)Uj!ix). ( 12 ) 

In order to remove potential bias from the vortex-only 

configurations being in the gauge where all elements are 
diagonal, we perform a random gauge transformation on 
them before smoothing. 


III. SMOOTHING 

Motivated by the recent results of Ref. |2B], showing 
that vortex-only configurations are able to reproduce dy¬ 
namical mass generation after 10 sweeps of cooling, we 
perform smoothing on our ensembles. We seek to quan¬ 
tify the extent to which the underlying vacuum structure 
of vortex-only and untouched configurations is similar af¬ 
ter smoothing. 

Smoothing is known to reveal a gauge field background 
consisting of a “liquid” of (anti-)instanton-like objects 
[351 136j . While instanton/anti-instanton pairs continue 
to annihilate under smoothing, isolated instanton like 
objects remain stable under further smoothing, as long 
as over-improved or highly improved actions are used. 
(Anti-)Instantons are associated with approximate zero- 
modes of the Dirac operator and thus through the 
Banks-Casher relation [28] contribute to a non-zero quark 
condensate, an order parameter for the dynamical break¬ 
ing of chiral symmetry. 

We will thus focus in particular on the instanton de¬ 
grees of freedom of our ensembles. We will smooth 
in two ways; cooling [37H1I] using an C>(a^)-three- 
loop improved action [3S1I1S] and over-improved stout- 
link (OISL) smearing [32I13 ]j an improved form of stout- 
link smearing tuned to retain instanton-like objects. 
Topological charge density is calculated using an 0(a^)- 
five-loop improved definition of the field-strength tensor 

m- 

In Fig. [^ we have plotted the ensemble-averaged ac¬ 
tion on our three ensembles as a function of smoothing 
sweeps. The untouched ensemble behaves as expected; 
the first few sweeps of smoothing rapidly removes action, 
as quantum fluctuations away from the classical solution 
are removed. Then, after about 10 sweeps of cooling or 30 
of OISL smearing, a stable background of instanton-like 
objects is revealed, and action is very slowly removed 
from the lattice as these can only disappear via pair- 
annihilation. 

The vortex-removed ensemble shows a very rapid loss 
of action, well below the untouched, and settles at a much 
lower value at high amounts of both cooling and OISL 
smearing. This corresponds to observations in Ref. [55] 
that vortex removal destabilises instanton-like objects, 
and so after smoothing we have produced a mostly trivial, 
low action background. 

The vortex-only ensemble, however, behaves markedly 
differently. The vortex-only configurations start very far 
from a classical solution to the QCD equations, and so 
have high action, decreasing only slowly under smooth¬ 
ing. Eventually, a rapid loss of action occurs. After 
around 20 sweeps of cooling, or 200 of OISL smearing, 
the rate of action loss decreases to be comparable to the 
untouched ensemble, and remains so up to high levels of 
smoothing. This corresponds to a dramatic change in the 
gauge field structure; it seems that this level of smooth¬ 
ing is required to produce a background of instanton- 
like objects on the vortex-only ensemble, which is then 

































4 



FIG. 3. A log plot of the ensemble-averaged action on untouched (squares), vortex-only (circles) and vortex-removed (crosses) 


ensembles under cooling (a) and over-improved stout-link smearing (b) 


very stable under smoothing, analogous to that found 
in the untouched case. Just as vortex removal desta¬ 
bilised instanton-like objects, the centre vortex structure 
alone contains the essential information for these objects, 
which are then revealed through the smoothing proce¬ 
dure. 

The evolution of the ensemble-averaged absolute value 
of the topological charge on the three ensembles is shown 
in Fig. 1^ Again, the untouched ensemble behaves as 
expected. Initially, the integrated topological charge is 
poorly defined on rough configurations, but after smooth¬ 
ing settles on a stable value. Once our configurations 
have been smoothed enough to resemble an instanton 
liquid, instanton-anti-instanton pairs disappear from the 
lattice only by pair annihilation, which has no net effect 
on the configuration topological charge, and so this value 
is stable up to very high levels of smoothing. 

Consistent with the action, the vortex-removed results 
show an almost empty background of very low topolog¬ 
ical charge after smoothing. The removal of centre vor¬ 
tices has destabilised otherwise topologically non-trivial 
objects, resulting in their subsequent destruction by the 
smoothing algorithm. The vortex-removed results do 
not, however, settle on no topological charge, indicat¬ 
ing imperfections in the vortex removal procedure have 
left some objects intact. 

In agreement with the action, the vortex-only results 
show a high degree of roughness until around 10 sweeps 
of cooling and 180 of OISL smearing, where they become 
much smoother, and display a similar result to the un¬ 
touched ensemble, settling on a stable value. 

To study the instanton degrees of freedom, we directly 
examine the gauge fields using the method of Ref. |12] . 
We search the lattice for local maxima of the action den¬ 


sity, around which we fit the classical instanton solution; 
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So{x) = C ^ 


TT^ ((x - Xo^ + P^)^’ 


(13) 


with p the instanton radius, and xq the location of the 
centre of the instanton, not restricted to lie on lattice 
sites. The scale parameter ^ is introduced to ensure we 
fit to the shape of the action density, not the magni¬ 
tude. We note that this method requires a relatively 
smooth background; UV noise will obscure the shape of 
instanton-like objects, and lead to fitting to a number 
of “false positives”; local maxima of the action originat¬ 
ing due to short-range noise, not larger-scale underlying 
vacuum structure. 

We have plotted the ensemble-average number of local 
maxima of the action density as a proxy for the number 
of instanton-like objects on the lattice as a function of 
smoothing in Fig. The rapid decrease after a small 
amount of smoothing in the untouched case corresponds 
to the appearance of a gauge-field background dominated 
by instanton-like objects; at small levels of smoothing 
local maxima are largely due to short-range noise. Af¬ 
ter this, the remaining objects are highly stable under 
smoothing. 

As expected, the vortex-removed ensemble becomes 
almost empty under smoothing; 2 orders of magnitude 
fewer objects are found under high levels of smoothing. 
Removing centre vortices from the lattice has correspond¬ 
ingly destabilised instantons, leading to their destruction 
by smoothing algorithms, although some remain. 

In the vortex-only case, the results seen for the inte¬ 
grated action and topological charge hold true for the 
action density; an initially noisy background, shown by a 
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FIG. 4. The ensemble-averaged absolute value of the topological charge on untouched (squares), vortex-only (circles) and 
vortex-removed (crosses) ensembles under cooling (a) and over-improved stout-link smearing, (b) 



FIG. 5. The ensemble-average number of local maxima of the action found on untouched (squares), vortex-only (circles) and 


vortex-removed (crosses) ensembles after cooling (a) and over-improved stout-link smearing, (b) 


high number of local maxima of the action density, comes 
to resemble the instanton configuration seen on the un¬ 
touched ensemble after a “turning point” of 10 sweeps of 
cooling and 180 of OISL smearing, thereafter remaining 
stable under smoothing. 

In Fig. 1^ we have plotted the ensemble-average fit¬ 
ted radius of instanton candidates. The untouched case 
shows the expected shape [H] ; an initial drop correspond¬ 
ing to the removal of false positives due to short range 
noise, followed by a small, steady increase as the gauge 
field background becomes dominated by instanton-like 
objects. The vortex-only result, again, shows an initially 


noisy background, evinced by the large number of small 
objects, transforming into an instanton background al¬ 
most identical to that in the untouched case by around 
40 sweeps of cooling or 200 of OISL smearing. The sim¬ 
ilar slopes beyond this points means a small shift in the 
x-coordinate would place this trajectory on top of the 
untouched trajectory. In the vortex-removed case, the 
average instanton radius becomes very large after the ini¬ 
tial cooling. This likely corresponds to an almost empty 
background, where the instanton fitting algorithm has at¬ 
tempted to fit over large distances with an almost smooth 
gauge field. After a large amount of cooling, and the large 
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(a) (b) 


FIG. 6. The ensemble-average radius (p) of instanton candidates found on untouched (squares), vortex-only (circles) and 
vortex-removed (crosses) ensembles after cooling (a) and over-improved stout-link smearing, (b) 


objects have been smoothed away, the vortex-removed 
ensemble shows a small, stable number of small instan- 
tons, most likely due to Gribov copy induced inefficien¬ 
cies in the vortex removal procedure failing to eliminate 
all centre vortices. 

To further elucidate the situation, we compare the fit¬ 
ted radius (p) and measured topological charge at centre 
{q{xo)) of instanton candidates to the theoretical rela¬ 
tionship, 

q{xo) = Q (14) 

where Q = =Fl for an (anti-)instanton. This is plotted for 
a typical configuration in Fig. [7] after 10 and 40 sweeps 
of cooling. Other configurations examined show similar 
behaviour. All instantons found on the untouched con¬ 
figuration very closely follow the theoretical relationship 
after 10 sweeps of cooling, and continue to do so at 40, 
although their number has been greatly reduced by pair 
annihilation. While the vortex-only configuration after 
10 sweeps of cooling has a large number of instantons ly¬ 
ing close to the theoretical relationship, there are also a 
large number of small objects lying far from the relation¬ 
ship. After 40 sweeps, only instanton-like objects remain. 
It seems that while instanton-like objects are present on 
the vortex-only configurations after 10 sweeps, it takes 
far longer than in the untouched case for smoothing to 
remove noise and thus create an instanton dominated 
background. Notably, there is no direct correlation be¬ 
tween the instanton-like objects found in the untouched 
and vortex-only cases. While we have recreated a similar 
instanton-liquid like background in the vortex-only case, 
we have not recreated specific individual objects. By con¬ 
trast, the vortex-removed background is almost empty 
after 40 sweeps, with only a small number of objects 
remaining. The way the smoothing algorithm destroys 


these otherwise non-trivial objects can be seen at 10 
sweeps; these objects are enlarged, resulting in a smaller 
topological charge at centre, until being smoothed away. 
This explains the observation of Fig. that at intermedi¬ 
ate ranges the objects on vortex-removed configurations 
are very large, before the average size shrinks dramati¬ 
cally. 

We have illustrated the qualitatively similar vortex- 
only and untouched backgrounds through the topological 
charge density in Fig.|^ We have plotted a level set of the 
topological charge density after 40 sweeps of cooling on 
a single configuration, both in its untouched and vortex 
only states. In both cases instanton-like objects are ap¬ 
parent, appearing as regions of high topological charge 
density. There is no 1-1 correspondence, however, be¬ 
tween the objects found on the untouched configuration 
and on the vortex-only one. 

The mechanism by which dynamical chiral symmetry 
breaking is apparent on vortex-only configurations after 
10 sweeps of cooling is thus made clear; although there 
is still a large amount of noise present, instanton-like ob¬ 
jects have begun to appear. Finding instanton-like ob¬ 
jects through direct examination of the gauge field re¬ 
quires a largely smooth background, while their presence 
can be felt throngh dynamical chiral symmetry breaking 
manifest in the overlap qnark propagator much earlier. 

Combined, the topological charge, action, and instan¬ 
ton characteristics paint a unified picture; under vor¬ 
tex removal, otherwise topologically non-trivial objects 
are destabilised, resulting in a trivial background un¬ 
der smoothing. By contrast, the vortex-only results, de¬ 
spite consisting solely of the centre degree of freedom of 
the original ensembles, are able to create the same long- 
range gauge field characteristics as seen in the untouched 
case after smoothing. The centre vortices, it seems, con¬ 
tain the “seeds” of instanton-like objects - the informa- 
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(a) 


(b) 




(c) (d) 




(e) (f) 


FIG. 7. The values of the instanton radius, p, against the topological charge at the centre, q{xQ), on a representative conhgu- 
ration. Results are compared to the theoretical relationship between the instanton radius and topological charge at the centre, 
Eq. (141 (solid lines), and the dislocation threshold, 1.97o (dashed line). Results are shown on a typical configuration for the 
untouched (top), vortex-only (middle), and vortex-removed (bottom) cases at 10 (left) and 40 (right) sweeps of cooling. Note 
the larger scale for the vortex-removed case after 40 sweeps of cooling. 


























FIG. 8. A visualisation of the topological charge density after 40 sweeps of cooling on a representative configuration, in both 
untouched (a) and vortex-only |(b)| form. Contours of positive topological charge density are plotted in yellow, and negative in 
blue, for a level set. 


tion necessary to recreate the topological structure of 
the gauge field background. Notably, both methods of 
smoothing have produced similar behaviour; the centre 
vortex degree of freedom robustly contains all informa¬ 
tion about the long-range vacuum structure, regardless of 
the smoothing method used to expose it. The measures 
of vacuum structure we have considered are also in agree¬ 
ment as to the amount of smoothing required for this to 
happen; it seems that at around 10 sweeps of cooling 
or 180 of OISL smearing, instanton-like objects have be¬ 
gun to appear on vortex-only backgrounds, and that by 
around 40 sweeps of cooling or 200 of OISL smearing the 
resulting vacuum structure closely resembles that of the 
untouched on a macroscopic level, after the same amount 
of smoothing. 

Having established the similarity of results under cool¬ 
ing and OISL, we will henceforth restrict calculations to 
cooled ensembles. 


IV. STATIC QUARK POTENTIAL 

Having seen the creation of a gauge field background 
on vortex-only configurations similar to the untouched 
case after smoothing, we now turn our attention to the 
ability of this background to recreate confinement. It 
has long been known [521 SS] that while vortex-removal 
results in the removal of the string tension, and thus con¬ 
finement, vortex-only ensembles, without smoothing, can 
reproduce only approximately 66% of the string tension. 

We thus examine the static quark potential on un¬ 
touched and vortex-only ensembles under cooling, plot¬ 
ted in Fig. 1^ At long distances, we have fit to the form 
V (i?) = aR + c, with fit parameter a summarised in 
table |T[ At 10 sweeps of cooling, results at long range 
are consistent with the unsmoothed case; the vortex-only 
string tension is only 67% of the untouched. However, 
after 40 sweeps of cooling, this picture has changed; the 
vortex-only ensemble now reproduces 93% of the string 



FIG. 9. The static quark potential on untouched (blue) and 
vortex-only (green) configurations after 10, 40, and 80 sweeps 
of cooling. Note that the 80 sweeps untouched potential is 
hidden behind the 80 sweeps vortex-only potential. 


tension found in the untouched case. This is accompa¬ 
nied by a decrease in the string tension found in the un¬ 
touched case, to approximately half its initial value. By 
80 sweeps, the percentage of the untouched string ten¬ 
sion reproduced by the vortex-only ensembles has risen 
to 97%. It seems that after smoothing the degrees of 


TABLE I. Value of the string tension (lattice units) found on 
untouched (UT) and vortex-only (VO) ensembles after 10, 40, 
and 80 sweeps of cooling. 

Cooling Sweeps crux crvo 

lO 1.120(4) 0.751(5) 

40 0.641(2) 0.592(2) 

80 0.423(1) 0.409(1) 
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freedom responsible for heavy quark confinement are also 
identical on the untouched and vortex-only ensembles. 

We note that while 10 sweeps of cooling are sufficient 
to recreate dynamical chiral symmetry breaking on the 
vortex-only configurations, 40 are required to recreate 
confinement. This difference in scales suggests an in¬ 
creased sensitivity to roughness when examining confine¬ 
ment. 


V. LANDAU-GAUGE OVERLAP QUARK 
PROPAGATOR 


We now examine dynamical mass generation, and cor¬ 
responding dynamical chiral symmetry breaking, on our 
ensembles, via the Landau-gauge overlap quark propa¬ 
gator. The Landau-gauge quark propagator provides a 
clear signal of the presence or absence of dynamical mass 
generation; the infrared enhancement (or lack thereof) of 
the Dirac scalar part of the propagator. We decompose 
the Landau-gauge quark propagator in momentum space 
into Dirac scalar and vector components as Hi ST] 


S{p) 


Z{p) 

H +M{p) ’ 


(15) 


where q is the kinematic lattice momentum |46j . M (p) is 
the non-perturbative mass function, and Z(p) contains 
the renormalisation information. We use the overlap 
fermion action, as it satishes the Ginsparg-Wilson rela¬ 
tion [?5], and thus provides a lattice-realisation of chiral 
symmetry. This superior sensitivity to chiral effects is 
key to examining the effects of gauge field topology. 

We use the fat-link irrelevant clover (FLIC) fermion 
operator mi EMI] as the overlap kernel with 

regulator parameter iriw = 1. We have considered 
three values of the overlap mass parameter, p = 0.004, 
p = 0.012, and p = 0.022, corresponding to physical bare 
quark masses of 12 MeV, 40 MeV, and 70 MeV respec¬ 
tively. We rotate to Landau gauge using a Fourier trans¬ 
form accelerated algorithm [52], fixing to the 0{a?) im¬ 
proved gauge-fixing functional |53j . A cylinder cut [54j is 
performed on propagator data, and Z{p) is renormalised 
to be 1 at the highest momentum considered, p ~ 5.2 
GeV. 

Results are presented in Figs.j^ and[^ As expected, 
the mass function in the untouched case displays strong 
infrared enhancement; dynamical mass generation via 
dynamical chiral symmetry breaking. The cooling proce¬ 
dure, which has left the instanton structure of the con¬ 
figurations intact, has correspondingly retained dynami¬ 
cal mass generation on cooled configurations. Although 
there is a small decrease in the amount of dynamically 
generated mass from 40 to 80 sweeps, due to pair anni¬ 
hilation of instantons, the qualitative long-range features 
of the quark propagator remain intact. At short range, 
cooling has removed all UV noise, and so consequently 
the mass function sits on the input bare mass at large 
momenta. 


The renormalisation function on the untouched ensem¬ 
ble likewise displays the characteristic shape; dipping 
slightly in the IR, reaching a plateau in the UV. These 
observations remain true at all bare quark masses con¬ 
sidered. 

As found in Ref. |2^, the vortex-only results are able 
to reproduce both the mass and renormalisation func¬ 
tions of the untouched ensembles. The background of 
instanton-like objects found after smoothing on vortex- 
only configurations has successfully reproduced the fea¬ 
tures of the quark propagator. This concords with the 
results of Ref. m. which showed the overlap operator 
displaying a similar gauge field background to that found 
after smoothing. Although it requires a higher level of 
smoothing for topological objects to become apparent in 
the gauge field structure itself, the overlap is sensitive to 
their presence after just 10 sweeps. We have shown the 
equivalence of vortex-only and untouched results persist¬ 
ing up to a high level of smoothing, and at high bare 
quark masses. The similar background of instanton- 
like objects on vortex-only and untouched ensembles ob¬ 
served by directly examining the gauge helds equally pro¬ 
duces dynamical mass generation on both untouched and 
vortex-only configurations. 


VI. CONGLUSIONS 

We have examined the gauge held structure of conhgu- 
rations consisting solely of centre vortices under smooth¬ 
ing, comparing to the original conhgurations under the 
same smoothing. Through study of the action, integrated 
topological charge, and instanton content of these con¬ 
hgurations, we have found remarkably similar structures 
emerging under smoothing, regardless of whether cooling 
or over-improved stout-link smearing is used. In both 
cases, after smoothing we obtain a gauge held structure 
approximating a background of instantons, stable under 
further smoothing. Under cooling, instanton-like objects 
begin to appear on the vortex-only conhgurations after 10 
sweeps, and after 40 sweeps the gauge held background 
of both the vortex-only and untouched ensembles consist 
of a dilute instanton liquid. 

This provides a mechanism for the hndings of Ref. [55] , 
that after 10 sweeps of cooling vortex-only conhgurations 
support dynamical mass generation, and thus dynamical 
chiral symmetry breaking. The information contained 
in the centre vortex degree of freedom is sufficient to 
recreate the underlying long-range structure of the QCD 
vacuum, revealed through smoothing. We have also con- 
hrmed that the gauge held background of vortex-removed 
conhgurations is unstable under smoothing. 

Upon examining the string tension, we have found that 
vortex-only conhgurations are, again, capable of repro¬ 
ducing the untouched results under smoothing. The sim¬ 
ilar long-range gauge held structure on vortex-only and 
untouched conhgurations after cooling has produced a 
similar conhning structure, as well as similar dynamical 
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FIG. 10. The mass (left) and renormalization (right) fnnctions on nntonched (bine squares), and vortex-only (green circles) 
ensembles after 40 sweeps of cooling, at physical bare qnark masses of 12 MeV &[g, 40 MeV @&|3, and 70 MeV 
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FIG. 11. The mass and renormalization functions on untouched (blue squares), vortex-only (green circles) ensembles after 80 
sweeps of cooling, at physical bare quark masses of 12 MeV 0 &[g, 40 Mev 0&I3, and 70 MeV 0&|§. 
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mass generation. 

Although we have found similar gauge field background 
structure on the untouched and vortex-only configura¬ 
tions after smearing, this naturally has removed the 
short-range information of the gauge fields. It would be 
desirable to have a method of smoothing that would allow 
us to gain smooth centre-vortex configurations without 
removing relevant short-range structure, allowing direct 
comparison to unsmeared configurations. It would be in¬ 
teresting to explore an SU(3) extension to the technique 
proposed for the SU(2) case in Ref. [55] . 

We have demonstrated how gauge-field smoothing can 
restore agreement between the characteristic features of 
confinement and dynamical chiral symmetry breaking ob¬ 
served on untouched and vortex-only ensembles. There 
are two possible interpretations of this result. In the first 
scenario we can consider that the discrepancies are asso¬ 
ciated with gauge-field roughness which is removed under 
smearing or cooling, thereby restoring agreement. An al¬ 
ternative interpretation is that the essential fundamental 
feature of the QCD vacuum is the thick centre vortex, 
rather than the thin centre vortices that are identified on 
the lattice. Then, smoothing would be primarily under¬ 
stood as a process of creating thick centre vortices from 
the identified thin centre vortices. The characteristic size 
of thick vortices is around 1 fm; we have found that 10 
sweeps of cooling, corresponding to a cooling radius of 
1.2 fm, is sufficient to recreate dynamical chiral symme¬ 
try breaking from a centre vortex background. It is of 
note that while vortex-only configurations consist only 
of thin centre vortices, the vortex removal procedure is 
sufficient to remove thick centre vortices from the lattice. 
This would explain the result that while vortex removal 
is sufficient to spoil dynamical chiral symmetry breaking 
and confinement, smearing is required to reproduce these 
from vortex-only configurations. 

The two interpretations we have postulated can be dis¬ 
tinguished by their behaviour under the approach to the 
continuum limit. Under the first scenario, to ameliorate 
the roughness of the vortex-only fields it would be suffi¬ 
cient to hold the number of smoothing sweeps fixed as the 
lattice spacing tends to zero. Under the second scenario, 
it would be necessary to increase the number of smooth¬ 
ing sweeps used as the lattice spacing decreases, such 
that the physical smearing distance is fixed to preserve 
the scale corresponding to the characteristic thickness of 
centre vortices. This will be the subject of future work. 
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Appendix: MCG fixing coefficients 

For the Maximal Centre Gauge fixing procedure, we 
wish to maximize the quantity 

i?iocai(a^') = EjTrU(x')C/^(x')P+ (A.l) 

at a given lattice site x', for a local gauge transformation 

U(x'). 

This is achieved by restricting U(x') to be in one of 
the SU(2) subgroups of SU(3). Then, we parametrize an 
SU(2) matrix [U(x')]su( 2 ) by 

[U(x )]gu( 2 ) ~ ^4 I i^i fJj, (A.2) 

and embed the resulting SU(2) matrix in one of the 3 
SU(2) subgroups of SU(3). This allows us to rewrite the 
gauge fixing condition in terms of the SU(2) parametriza- 
tion as 

4 ^ 4 

-^local(^ ) ~ ^ ^ ^ ^ ('^■^) 

ij = l i 

Where are elements of a real, symmetric matrix, bi a 
real vector and c a real constant, all dependent only on 
U^{x') and U^{x’ — jl). Note that we suppress the sum 
over directions, /i. Once the coefficients are defined the 
local quantity i?iocai(a^0 be maximised according to 
the method of, e.g., Ref. (33). 

The gauge fixing coefficients a^, bi, and c are defined 
in terms of the gauge links U^{x') and U^(x' — /t), for 
any given SU(2) subgroup of SU(3). 

It is convenient to define 

U = U^{x '), 

V = Uf,{x - fi ), 

(A.4) 

Then we define complex 2x2 matrices, U and V, as the 
components of U and V corresponding to the given SU(2) 
subgroup, and complex scalars, u and v, as the remaining 
diagonal element. The other elements of U and V do not 
contribute to the coefficients, and so for convenience are 
set to 0. 

As an example, for the first SU(2) subgroup, we set 
u=( ^ 0 y (A.5) 

\ 0 0 u / 


and 
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V = 


0 0 V ) 


(A.6) 


Once we have defined U, u, V, and v, the values of the 
gauge-fixing coefficients depend only on these, and so the 
expressions below are valid for all SU(2) subgroups. 


Then we have for the real symmetric matrix a^, 


aii= |Ui2|' + |U2i|" + Ui2U;i -bUaiUts + |Vi2p + IV 21 I" + Vi2V;i + V2iVt2 , (A.7) 

012 = 2* (Ui2U*i - U2iUt2 + Vi2V*i - V2iVt2) , (A.8) 

ai3 = +UnU*,2 + Ui2Uti + UnU^i + U2iUJi - V 22 IJI 2 - Ui2U^2 - U22U^i - U2iU;2 

+ VnVt2 + Vi2Vti + VnV*i + V2lVti - V22Vt2 - Vi2V*2 - V22V*1 - V2 iV;2 , (A. 9 ) 

014 =i(UnUj 2 - Ui 2 Uti + UiiU*i - lJ 2 llJn + U 22 Ut 2 - Ui 2 U ;2 + U 22 U *1 - U 2 iU *2 

- ViiVt2 + Vi2Vti - ViiV*i + V2iVti - V22Vt2 + Vi 2V2*2 - V22V*i + V2iV;2 ) , (A-IO) 

022 = |Ui2p + |U2l|"-Ui2U;i - U2lUt2 + | V 12 | ^ + | V 21 ^ - V 12 V*i - V2lVt2 , (A.ll) 

023 =i(-UiiUl2 + UiaUti + UiiU^i - U21UI1 + U22Ut2 - Ui 2 U ;2 - U22U;i + U21U22 


-ViiVt2 + Vi2Vti + ViiV*i - V2lVti + V22Vt2 - Vi 2V*2 " V22V*i + V 21 V*2 ), (A.12) 

024 = +UiiU *2 + Ui 2 Uti - UiiU^i - U 2 lUti + U 22 Ut 2 + Ui 2 U ;2 " U 22 U;i - U 2 lU ;2 

- ViiVt2 - Vi2Vti + ViiV^i + V2iVti - V22Vt2 - Vi2v; 2 + V22v;i + V2 iV;2 , (A.13) 


033= |Uii |2 + |U22|"-UiiU; 2 - U22Uti + |Vii |2 + |V22|"-ViiV; 2 - V22Vti, (A. 14 ) 

034 = 2 * (- UiiU ^2 + U 22 Uti + Vii V *2 - V 22 Vti ) , (A. 15 ) 

044= |UiiP + |U22p+UiiU^2 +U22Uti + |Vii |2 + |V22p+ViiV22 + V22Vti , (A. 16 ) 

and for the real vector 6 ^, 

bi=i (-U12U* + U*2U - U21U* + U^iu + V12V* - Vt2V + V21V* - V^iv) , (A. 17 ) 

62 = U12 u* + Ut2 u - U21 u* - u - V12 v* - V:;;2 V -f V21V* -h v, (A. 18 ) 

b 3 = i (-U11U* + Ujiu + U22U* - U; 2 U + Viiv* - Vtiv - V22V* + V* 2 v) , (A. 19 ) 

&4 = Uii U* + UJi u + U22 U* + U^2 u + Vii v* + VJi V + V22 V* + V*2 V, (A. 20 ) 

and 

c=|up + |vp. (A. 21 ) 
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